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I show that for a fi>type organic (BEDT-TTF)2X molecular crystal, a superconducting state with 
T c ~ 10 K and gap nodes on the Fermi surface can be caused by short-ranged antiferromagnetic spin 
fluctuations. Using a two-band description for the anti-bonding orbitals on a BEDT-TTF dimer of 
the K-type salt, and an intermediate local Coulomb repulsion between two holes on one dimer, the 
magnetic interaction and the superconducting gap-function are determined self consistently within 
the fluctuation exchange approximation. The pairing interaction is predominantly caused by inter- 
band coupling and additionally affected by spin excitations of the quasi one-dimensional band. 



Electronic mechanisms for superconductivity are of 
conceptual interest on their own and allow, in princi- 
ple, large superconducting (SC) transition temperatures 
since electronic energy scales easily exceed phononic en- 
ergies. In addition to the copper-oxide based high- 
temperature superconductors, heavy fermion supercon- 
ductors like UPt 3 and UBe 13 , and the ruthenate com- 
pound Sr 2 Ru04, organic superconductors are candidates 
for a pairing mechanism which originates in strong elec- 
tronic correlations in the normal state. 

The organic molecular crystals k-(BEDT-TTF) 2 X 
(in the following abbreviated as k-(ET) 2 X), are very 
anisotropic quasi two-dimensional superconductors with 
transition temperatures up to around 10 K, depending on 
pressure and the ion X, which can, for example, be I3, 
Cu[N(CN) 2 ]Br or Cu(SCN) 2 §. Shubnikov - de Haas 
experiments have established the existence of a well de- 
fined Fermi surface in these materials m, demonstrat- 
ing the Fermi liquid character of the low energy quasi- 
particles. Below T c , the low temperature 13 C-NMR spin- 
lattice relaxation rate varies as T 3 the electronic 
specific heat coefficient exhibits a magnetic field depen- 
dence 7 oc f? 1 / 2 |?J and the thermal conductivity K is 
linear in T ||; evidences for nodes of the SC gap func- 
tion on the Fermi surface. All these results suggest that 
superconductivity is caused by a highly anisotropic and 
likely non-phononic pairing interaction. An electronic 
mechanism is also supported by the absence of an iso- 
tope shift of T c due to the C=C and C-S phonon modes, 
as observed by 34 S and 13 C= 13 C isotope substitution in 
the central part of the ET-molecule j^] . 

It has been pointed out earlier that many anomalous 
normal state properties of layered organics strongly sug- 
gest the important role of electronic correlations 10-1^] 
and that some, but not all, of these anomalies are in 
striking similarity to cuprate superconductors p3|. Par- 
ticularly magneto-transport, the thermoelectric power, 
and the uniform magnetic susceptibility display similar 
behavior. Another resemblance is their proximity to an 
antiferromagnetic (AF) state, even though, in organics, 



frustration and strong in-plane anisotropics demonstrate 
differences between both systems. Particularly, due to 
the, by now established, dominant d x 2_ y 2 order parame- 
ter component O], high-temperature superconductivity 
in cuprates is believed to be dominated or entirely caused 
by an electronic mechanism. The spin fluctuation sce- 
nario [l5| , [T^] , which predicted the d x 2_ y 2 state, is a very 
promising approach, because it also offers an explana- 
tion for a variety of anomalous normal state properties of 
cuprates. The unfrustrated AF nearest neighbor coupling 
and the specific shape of the Fermi surface are important 
prerequisites for a spin-fluctuation induced SC state with 
d x 2_ y 2 symmetry in cuprates. In contrast, in organics, 
the magnetic coupling is frustrated, due to the underly- 
ing anisotropic triangular lattice | fl3|| . In addition, their 
Fermi surface consists of two disconnected pieces 
Therefore, it is not at all obvious whether, despite these 
numerous similarities to cuprates, spin fluctuations can 
bring about superconductivity in the organics. 

In this paper I show that spin fluctuations can cause 
superconductivity with T c sa 10 — 15K in K-type or- 
ganic molecular crystals fT^ ] and discuss similarities be- 
tween the inter-band coupling and the hot spot scenario of 
cuprates |jl5|,[l8). Gap nodes on the Fermi surface, accom- 
panied by changes in sign of the order parameter, result 
because of zone boundary inter-band pairing combined 
with a intra-band coupling of a quasi one-dimensional 
(quasi 1-D) band. This is in agreement with NMR and 
specific heat measurements and should be observable in a 
corresponding phase shift experiment. Furthermore, the 
strong inter-band pairing interaction causes an opposite 
sign of the gap functions of the two bands. 

Within a given conducting layer, the unit cell of k- 
(ET) 2 X is occupied with six electrons in the outer elec- 
tronic shells and consists of four ET- molecules, which 
can be considered as single electronic degrees of free- 
dom [ fi9| . The four molecules are arranged in two dimers 
(see Fig. Qa). Due to the large intra-dimer transfer in- 
tegral, t\ pa 250 meV, the two bonding bands are com- 
pletely occupied and shifted to binding energies ~ 2t\ 
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below the Fermi energy. Only the two remaining anti- 
bonding bands cross the Fermi surface; those are half 
filled with the remaining two electrons. 
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FIG. 1. (a) Spatial arrangement of the ET molecules within 
the planes. The two dimer types are indicated by the dif- 
ferent tilting of the molecule pairs, which are, after dimer- 
ization, considered as a single electronic degree of freedom. 
The dashed (dotted) lines indicate the unit cell for tz 7^ ty 
(*3 = *3')- 00 Fermi surface for the case £3 7^ ty , consist- 
ing of a quasi 1-D band and a hole pocket closed around the 
Z-point. (c) Fermi surface for the case £3 = ty. Note that the 
hot spots due to nearest neighbor AF correlations correspond 
to the gaps between both Fermi surface parts in (b). 

In what follows, I consider only the two anti-bonding 
bands near the Fermi surface (within a hole picture) and 
neglect the bonding bands. Due to the strong intra- 
molecular Coulomb repulsion (Uet > 1 eV), a doubly oc- 
cupied dimer will likely distribute the two holes on both 
molecules of the dimer, causing an effective correlation 
energy U ~ 2t± determined by the bonding-anti-bonding 
band splitting O] . This leads to the following two band 
Hubbard Hamiltonian 
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Here, ct is the creation operator of a hole with spin 
a in the anti-bonding band of the Z-th dimer (J, £ 
{1,2}) in the i-th unit cell and nn a = <\i„CilcT- Uj 
and tij are inter-dimer hopping elements between dimers 
of the same and of different type, respectively. Their 
Fourier transforms *k and <k with k = (k y ,k z ) de- 
termine the band-structure of the two bands £k,± = 
tk ± tk- Various different tight binding parameteriza- 
tions have been proposed for *k and PJl9|-pl[. I 
will use the dispersion relation *k = 2t2COs(k v ) and 



? k = 2cos(k z /2)y/tf+q, + 2*3*3/ cos(k y ) [|J2|], which 
can be derived from a four band model up to leading or- 
der in t-2.jt\. Due to the non-symmetrical position of the 



dimers and the ET-molecule structure itself, the three 
hopping elements are slightly different. The Fermi sur- 
face, for t% — 45 meV, £3 = 60 meV and ty = 65 meV, 
shown in Fig. [l]b, reproduces de Haas - van Alphen mea- 
surements |j| and electronic structure calculations ]l9| . 
Note, the two pieces of the Fermi surface are disconnected 
since £3 7^ ty . 

For the investigation of an electronic pairing state, one 
has to determine, the momentum and frequency depen- 
dence of the pairing interaction which is responsible for 
the magnitude, A^, |, of the coordinate space compo- 
nents of the gap function All oc {cmCji'l) an d the rela- 
tive phase relations the Am establish. The latter effect is 
determined by the maximum gain of condensation energy 
on the Fermi surface. In order to have a quantitative ac- 
count for this interplay of magnetic correlations, Fermi 
surface shape and inter- and intra-band coupling, I use a 
self-consistent summation of bubble and ladder diagrams 
(fluctuation-exchange approximation [^,^3|) within the 
Nambu-Gorkov description of the SC state. Both the 
pairing interaction and the SC gap function will be de- 
termined self consistently. Whether superconductivity 
occurs and whether the gap function possesses nodes on 
the Fermi surface will therefore be the results of our cal- 
culation. 

The fluctuation-exchange approximation has been suc- 
cessfully used for the investigation of superconductivity 
in one band models relevant for high temperature super- 
conductors p2]-p^|. For completeness, I give the set of 
equations for the two band Hamiltonian, Eq. |l|, within 
the SC state. The effective interactions of the system, 
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±^-[(l + U X c (q))- 1 -l]x C (q), (2) 



are (2x2) matrices in the dimer representation with 

particle-hole bubble Xu^ (?) = - J2k (k)Gi>i(k + q) + 
{—)Fw (k)Fi'i(k+q) for spin- and charge- type excitations. 
Here Gw (k) and Fu> (fc) are the normal and anomalous 
Green's functions with q = (q,iz/ n ), k = (k, iu) n ) as 
well as J2 k = T/NJ2k n - ^> *l are the two-dimensional 
wave vectors and the Matsubara frequencies are given 
by uj n = (2n + 1)ttT for fermions and v n = 2nirT for 
bosons. T is the temperature and N the number of unit 
cells. From the effective interaction, Eq. 0, in the dimer 
representation, one can obtain the intra and inter-band 
interactions V xx ,{q) with A, A' G {+, — }. The self energy 
in Nambu representation is assumed to be diagonal in 
the band representation: 



t x {k) = Y x (k)r a + X x (k)r 3 + 1> x (k)i 



(3) 



Here, Y x (k) = iu> n (l — Z x (k)) and X x (k) are the diag- 
onal and $\(k) the off -diagonal components of the self 
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energy in Nambu representation, respectively, and can be 
determined from the Eliashberg equations, 



Y x (k) = J2 V+ X ,(k - k')iu n ,Z x ,{k')lD x ,{k') 
k',X' 

x x (k) = J2 v \y( k - k ') +Xx>{k'))/D x ,{k') 

$A 
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( fc ) = E [ F A> (* - *') + U] *v (k')/D x , (k') , (4) 

it', A' 



with denominator D x (k) = (iuj n Z x (k)) 2 — (ek,A + 
X x (k)) 2 — Q x (k) 2 . The off-diagonal self energy, which 
signals superconductivity, determines the gap function 
A^(fc) = &\(k)/Z x (k). After analytical continuation to 
the real axis,zw„ — » uj + iO + 1 this set of coupled equations 
is solved self consistently using the numerical framework 
described in Ref. I25|l. 



of the two bands: Ak.+ ~ — Ak.-. Furthermore, the pro- 
nounced anisotropy of the pairing interaction, peaked at 
q = (±7r, 0), will cause a substantial momentum depen- 
dence of Ak,±. 
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FIG. 2. Effective pairing interaction for intra and in- 
ter-band excitations along the high symmetry lines of the 
Brillouin zone. Note the pronounced peak at q = (n, 0) for 
inter-band excitations. 

In Fig. 0a, we show the unit cell and in Fig. [TJd the 
corresponding Fermi surface, determined by £k,±- The 
Fermi surface segments of the two bands appear as one 
quasi 1-D and one hole pocket part, as seen in magneto- 
oscillation experiments |^,^|. This Fermi surface gen- 
erates 2kp intra-band transitions within the quasi 1-D 
band. In a real space picture, this corresponds to a cou- 
pling between dimers of the same type. In addition we 
expect strong couplings J3.3/ oc t 2 3 ,/U between dimers 
of different type. The numerical analysis of the above set 
of equations shows that this latter effect, which causes 
a magnetic inter-band coupling, dominates. Correspond- 
ingly, the inter-band effective interactions V^rj__(w) are 
large compared to V^ ++ (lu) — V^__(lu) and peaked at 
q = (±7r,0). This result is shown in Fig. 2, where we 
plot V~ xx , (u> = 0) along the high symmetry lines of the 
Brillouin zone (BZ). Due to this strong inter-band cou- 
pling we expect a relative sign —1 for the gap functions 



FIG. 3. Polar plot of the superconducting gap functions 
Ak,± of the two bands along the Fermi surface. For a given 
angle tan _1 (fc z /k y ) with Fermi vector (k y ,k z ), the amplitude 
of the gap corresponds to the distance from the origin. Gap 
nodes occur along the k y - and fc z -axis. The suppressed gap 
close to the diagonal is due to the absence of a Fermi surface. 
If tz = ty this region vanishes with equal gap amplitude at 
the meeting point of the two Fermi surface parts. 



Our numerical results for the anomalous self energy 
$k,±(w = 0) oc Ak.± are shown in Fig. 3. We find 
a stable SC solution, caused by AF spin fluctuations 
at a temperature T — 8 K which exhibits zeroth of 
the gap-function along the lines k y , z = 0. The lead- 
ing harmonics of the gap-function are given by Ak,± ~ 
A^sin(fc !; /2)sin(fc z /2) with A^ = -0.9 meV and A = 
1.5 meV. For T = 15 K, within our numerical accuracy, 
the gap function vanishes, consistent with an estimate for 
the superconducting transition temperature of the order 
of 10 — 15 K, in fair agreement with the experiment. As 
shown in Fig. lc, this gives gap nodes for both pieces 
of the Fermi surface. The symmetry of this supercon- 
ducting state is d xy which, in addition, is out of phase 
between the two parts of the Fermi surface. 

The physical origin of this pairing state and the close 
resemblance to the d x 2_ y 2 state in cuprates can be un- 
derstood if, for illustrating purposes, one makes the as- 
sumption t 3 — ts> for the hopping elements between dif- 
ferent dimer types. Now, the two dimer types, which 
differ in their orientation of the molecule pairs, arc indis- 
tinguishable and the unit cell reduces to that indicated 
by dotted lines in Fig. ^a. Assuming furthermore that 
the dimers are arranged on a square lattice, the BZ is 
doubled and rotated by tt/4. The resulting Fermi sur- 
face is shown in Fig. [lj:. One easily recognizes that the 
two branches of the Fermi surface of Fig. ^b correspond, 
after down- folding into the reduced BZ, to states sepa- 
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rated by the dashed line in Fig. [j]c. Since in coordinate 
space the dominating AF coupling is between nearest 
neighbors, it causes peaks in the magnetic susceptibil- 
ity around (±7T, ±7r) of the extended BZ. A quantitative 
calculation shows that the additional anisotropy due to 
the hopping clement £2 causes the peaks at ±(tt,tt) to 
differ from those at ±(7r, — it). Nevertheless, (±7r, ±ir) 
arc the dominating momentum transfers causing anoma- 
lies for hot quasiparticles, located on Fermi surface seg- 
ments close the corresponding magnetic zone boundary 
(dashed line in Fig. |I]c) |j^,|26|], which lead to the inter- 
band coupling in the original BZ of Fig , [ji b. On the other 
hand, as discussed in detail in Refs. Jl5|Jl^] , a momen- 
tum transfer (±7r, ±7r) causes a superconducting state 
with gap function Ak oc cos k y — cos k z . Rotation by 
7r/4 and down folding into the reduced BZ of Fig. [l]b 
yields Ak,± oc ± sin(fc a /2) sin(fc z /2) which covers the 
general behavior of the numerical results presented in 
Fig. 3. Thus, subject to the modifications brought about 
by the gap between the two branches of the Fermi surface 
and the additional anisotropy of the underlying lattice, 
the origin of the pairing state in cuprates and organics 
is very similar. Due to the important role played by hot 
momentum states for various normal state anomalies of 
cuprates jl8],^6|, we expect that they are responsible for 
normal state anomalies in organics as well. 

In conclusion, I have shown that spin fluctuations are 
a promising candidate for the pairing interaction of the 
superconducting state of organic molecular crystals. I 
find a transition at 10 — f 5 K to a superconducting state, 
in which the calculated occurrence of nodes on the Fermi 
surface is in agreement with various experimental obser- 
vations. Furthermore, a close similarity of the origin of 
this pairing symmetry to the one in cuprate supercon- 
ductors was demonstrated. Therefore it is tempting to 
conclude that the occurrence of two distinct classes of 
quasiparticles, hot and cold, depending on the strength 
of the effective interaction, shown to be vitally impor- 
tant for the normal state anomalies of cuprates, are es- 
sential for the corresponding unconventional behavior in 
organics as well. This would demonstrate that the spin 
fluctuation model provides a general scenario for systems 
characterized by antiferromagnetically short ranged and 
over-damped spin excitations, independent of the specific 
details causing these collective spin modes. Nevertheless, 
it should be noted that the self consistent weak-coupling 
approach, used in the present paper, does not cover 
strong coupling effects like the formation of a frustration- 
induced spin liquid state nor does it provides the desir- 
able quantitative description of the low-frequency spin 
dynamics as reached in cuprates p5[ ; rather it is only 
expected to give qualitatively correct information about 
the occurrence and symmetry of a SC state. 
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